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Abstract—In this paper we present an alternative 
approach to addressing the problem of scattering 
reduction for radar targets, which have recently been 
dealt with by using the Transformation Optics (TO) 
algorithm which calls for the use of Metamaterials 
(MTMs) that are inherently narrowband, dispersive and 
highly sensitive to polarization as well as to the incident 
angle. The present design utilizes realistic lossy materials 
that can be conveniently fabricated in the laboratory, and 
are wideband as well as relatively insensitive to 
polarization and incident angle of the incoming wave. A 
modified interpretation of the TO algorithm is presented 
and is employed to design of RCS-reducing absorbers for 
arbitrarily shaped targets, and not just for canonical 
shapes, e.g., cylinders, for which cloaks have been 
designed by using the TO. The paper also briefly 
examines the topic of performance enhancement of 
absorbers by using graphene materials and embedded 
Frequency Structure Surfaces (FSSs). 


We begin by presenting the design procedure for 
planar absorbers, and then examine how well those 
designs perform for general-shaped objects. Finally, we 
discuss how the planar design can be modified by 
tailoring the material parameters of the coating for 
specific object shapes. A number of test cases are 
included as examples to illustrate the application of the 
proposed design methodology, which is a modification of 
the classical TO paradigm. 


Index Terms—Transformation Optics (TO), cloaking, 
RCS reduction, RFI, microwave absorber, EMI. 


I. INTRODUCTION 


The objective of this paper is twofold. First of these to 
review the basic principles of the Transformation Optics 
(TO) approach, also known as Transformation 
Electromagnetics (TEM) algorithm, which has recently 
surfaced as one of the most innovative techniques for 


designing a wide variety of electromagnetic devices 
including cloaks to render objects invisible. Our second goal 
is to present an alternative approach for designing absorptive 
coatings for scattering reduction, which differs from that 
used to design TO-based cloaks, and avoids the problems of 
narrow bandwidth and sensitivity to polarization and incident 
angle associated with the TO-based cloaks that we realize by 
using Metamaterials (MTMs). 


The TO algorithm for designing cloaks is unique in that it 
provides a systematic approach, which is very innovative, as 
well as elegant, and is markedly different from the 
techniques that have been employed heretofore prior to the 
advent of the TO, to design RCS-reducing absorbers, for 
example. 


(a) (b) 
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Fig. 1. (a) Physical and (b) virtual domains used in the TO algorithm. 


Il. FUNDAMENTALS OF TO 


The principle upon which the TO is based has been 
enunciated in a number of papers, dating back to almost fifty 
years ago, that have examined the behavior of Maxwell’s 
equations in a generalized curvilinear coordinate system. To 
explain the basic principle of the concept, let us consider two 
objects belonging to physical and virtual domains, and 
shown in Figs.1 (a) and (b), respectively. The medium 
parameters surrounding these objects, namely (€), w,) and (ez, 
u2), are also shown in Fig.1. 


A number of authors, prominent among them Pendry [4], 
Leonhart [10], Hao [12, 42, 43] and Werner [11], as well as 
several others [13-22], have presented the relationship 
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between the medium parameters in the two domains 
(physical and virtual) when we transform the geometry of 
object #2 in the virtual domain into that of object#lin the 
physical domain via coordinate transformation, under the 
physical constraint that the electric and magnetic fields in the 
two domains remain “invariant” to the transformation 
between the two systems. The relationship can be explicitly 
stated as: 


OX, Ox, Ox, 
OX, Oy, OZ 
aA OL 1 
OX> OV> OZ 
OZ, OZ, 02, 
OX Yə OZ 
(la) 
= A "Ê -AT 
«1 “det (A) 
= _ A+ fly “AT 
“1 = “Get (A) 
(1b) 


where A represents the Jacobian matrix relating the two 
domains (x, y and z can be any arbitrary curvilinear 
coordinate), € and mw represent the permittivity and the 
permeability of the corresponding mediums, respectively. 
Eqn. (1) enables us to navigate between the two systems and 
relate their material parameters in a unique, systematic and 
rigorous way. 


Pendry[4], Smith [3] and a number of other workers [7-22] 
have leveraged the fact that the medium parameters can be 
related via (1), in order to lay the foundations of the TO 
algorithm for designing cloaks, which render a target 
invisible when covered by using materials whose parameters 
are dictated by the TO. To explore the basic principles of the 
TO, we return to Fig.l and define the following task for 
ourselves: Design the cloak (1.e., a cover or a coating) for the 
PEC object in Fig.1 (a) such that it is invisible to an arbitrary 
incident field that impinges upon it. Note that no restriction 
is being placed on the frequency, polarization or the angle of 
incidence of the illuminating field in connection with this 
task at this point. 


To solve the problem posed above, we begin with a 
coordinate transformation, which morphs object #1, which is 
located in the physical domain, and for which we are trying 
to design the cloak, into object#2 residing in the virtual 
domain. Although this transformation is not obvious when 
the geometries of both the objects are totally arbitrary, it 1s 
nonetheless doable, at least theoretically. The caveat, though, 
is that the procedure provides us no guarantees that the 
materials parameters dictated by the transformation can be 
realized physically, and/or that they can be fabricated in 
practice to achieve cloak designs which satisfy the realistic 
specifications, such as small thickness, wide bandwidth, 
polarization insensitivity, etc. 


The realizability issue, alluded to above, becomes even 
more critical when we attempt to design an “invisible” cloak, 
which is highly sought-after by cloak designers. We will now 
explain why this is the case with a simple example shown in 
Fig.2(a). 

Let us assume, for the sake of convenience, that the target 
we wish to cloak is a cylinder. Following the TO paradigm, 
this problem has been extensively studied by a number of 
authors [1-22] and they have derived the material parameters 
for the cloak by invoking the TO, which makes the cloaked 
object disappear entirely (become invisible). In Figs.2(b) and 
(c) we plot the material parameters, presented in [9], that are 
required to make the cylinder invisible. 


(a) 
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Fig. 2. (a) Schematics for a cloaked PEC cylinder with R, = 2R,, (b) 
material parameters for an all-angle, all-polarization cloak, (c) material 
parameters for a normal-incident, TE-polarization cloak from [9]. 


We observe several things from Fig.2. First, we see that 
the material parameters are anisotropic, and this in of itself 
can be problematic when we attempt to realize them in 
practice, because there is no systematic method available for 
synthesizing them. The second thing we observe about the 
required u and e values is that they vary over a wide range, 
tending to 0 in some regions and œ in others. Once again, 
this type of inhomogeneous behavior and a wide swing in the 
required material parameters as functions of the radial 
distance make it very difficult to realize them. In fact, we 
must resort to using Metamaterials (aka 
artificially-engineered materials) that are notoriously 
narrowband, dispersive and lossy where they attempt to 
realize the above type of material values. The above 
undesirable attributes that degrade the performance of the 
cloak and render it unsuitable for most applications that 
involve scattering reduction for real-world targets. What 
exacerbates the problem even more is the fact that the 
thickness of the cloak is comparable to the wavelength, 
rather than being a small fraction of the same. As we well 
know, a thin coating is desired in most applications, e.g., 
when designing stealth targets for the radar world. 


At this point we return to the TO paradigm for cloak 
designs and scrutinize it carefully to see if we can thresh out 
the root causes of the difficulties that we have just identified 
with the TO-based design, where upon introduction of 
Metamaterials for physical realization of the cloak structure. 
Note that polarization-insensitive Metamaterial designs for 
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TO-based cloaks will be extremely challenging and the usual 
simplification is to limit the polarization to a certain linear 
polarization. This introduces polarization sensitivity of the 
cloaks. Another example is that typical Metamaterials are 
designed using periodic structures and the unit cell details 
are optimized for a single angle of incidence (usually normal 
to the interface). This introduces sensitivity to the angle of 
incidence. 


To identify the problem areas with the TO-based cloak 
designs, we turn to an alternate derivation of the material 
parameters in the context of Transformation 
Electromagnetics. Rather than relying directly on the 
Jacobian of the transformation, which relates the Physical 
domain to the corresponding virtual domain, we turn to the 
integral forms of Maxwell’s equations appearing below: 


(2a) 


(2b) 
Initially, we consider the relatively simple case where the 
two domains are simply related by scaling, say by a factor 
‘y’, as has been done in the previous TO-based cloak designs 
[9].The method we propose to relate the material parameters 
in the two domains is very general and is applicable to the 
case where the two geometries have arbitrary shapes, and are 
not simply related to each other by a scale factor, as we will 
discuss later in Sec. III. 


2c 


2b [12a 
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Fig. 3. TO-based cloak schematics and corresponding materials in the (a) 
physical geometry and (b) virtual geometry. 


To derive the material properties of the physical domain, 
from the assumed parameters in the virtual domain (€o, uo for 
free-space in this example), we turn to Fig.3 and (2). The 
next step, we choose to impose the condition that the fields 
(E1, H1) in the physical domain be “identical” to those in the 
virtual domain, 1.e., (Ez, Hz). To facilitate the imposition of 
this condition, we now discretize the two domains, as shown 
in Fig.4, by setting up a mesh to discretize the regions-2 and 
-3in both domains. We take advantage of the circular 
symmetry of the geometries in the two domains, and of the 
fact that the geometry of the PEC cylinder, located in 
region-2 in the virtual domain is simply a scaled-down 
version of the one in the physical domain (region-1), and 
note that this transformation preserves the azimuthal 
symmetries of the two domains and, hence, the mesh size in 


the azimuthal direction remains unchanged when we 
navigate between the two domains. However, we follow a 
different strategy in the radial direction, along which we 
impose the following three conditions: 


(1) The number of radial cells be identical in region-2 of 
the domains, which spans the radial distance b<r<c in the 
physical domain and a<r<c in the virtual domain (note we 
have chosen c, the outer radius of region-2 to be identical in 
both domains); 


(11) The cell size in the outermost boundaries of region-2 1s 
to be identical in the two-domains; 


(iui) The cell sizes in region-3 in both domains be 
identically equal in the radial direction, not only between 
c <r<c, but also beyond, i.e., c<r<oo, meaning all the way to 
infinity. 






Region 2: Cloak Ma 





(a) (b) 


Fig. 4. Cloak problem with mesh schematics: (a) physical and (b) virtual 
domains. 


There are logical reasons for imposing the above 
constraints, as we will now explain, before we proceed to 
derive the relationships between the material parameter 
distributions in the two domains. First, the TO must be 
applied to entire regions external to the PEC target—cylinder 
in this case—when relating the material parameters in the 
two domains, viz., b<r<oo in the physical domain and 
a<r<oo in the virtual model. Specifically, we cannot truncate 
the regions as we transform from one domain to another, 
without introducing discontinuities in the fields, and thus 
violating the premise of TO under which we are operating. If 
we use condition (iii), in (2a and 2b), we immediately see 
that if the two fields as well as the mesh sizes would be 
identical in the two domains. It follows, then, that the 
material parameters must be exactly the same in regions-3 in 
the two domains. Since we have chosen the material 
parameters in the virtual domain in region-3 to be free-space, 
1.e., these parameters are € ) and fo, then region-3 in the 
physical domain must correspond to free-space as well, as 
would the parameters of the external region c <r<oo. This is 
a crucial point, and it implies that the cloak in the physical 
domain has but a finite thickness, spanning the region b<r< 
c. This is obviously necessary in order for the cloak design 
to be practical, since we cannot accept a cloak design whose 
thickness is infinite, as it would be if we did not impose the 
equal-mesh-size condition in the two domains in region-3 
and beyond, 1.e., for r > c. 
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To solve this problem we turn to the differential form of 
the Maxwell’s Equations given below: 


1 
Pal aS 
at’ 
(3a) 
1 
H-Al=+—e-E-AS 
At 


(3b) 
where Z S, and t represent the perimeter of a cell, 
surface area of a cell and time step in FDTD, respectively. 
Note that we are considering only the £, and H; components 
for this 2D geometry, namely a cylinder, for which z is along 
the axis of the cylinder, and € and u in (3) are the appropriate 
elements of the and l tensors. 


Let us now turn to the regions-2 in the two domains, 
namely physical and virtual. Recall that we have imposed the 
condition that the cell size at the outermost boundary of this 
region be identical in the two domains, which guarantees that 
the transition of the material parameters would be smooth as 
we transition from region-2 to region-3. (Recall region-3 and 
beyond is free-space in both domains.) So, all that remains 
for us to do now is to determine the material parameters of 
the cloak region in the physical domain, which spans from b 
< r< c, by invoking the condition that the fields (E4, H1) of 
the physical domain be identical to the fields(E,, Hz) in the 
virtual domain (which we have chosen to be free space), as 
well as condition (1) on the number of cells in the two 
domains associated with region-2, namely that this number 
be identical in the two domains. At this point, we impose an 
additional condition, without loss of generality, that in the 
radial direction the cell sizes in region-2 in the physical 
domain be all equal, as we go from r = b to r = c. We also 
choose the cell size in the radial direction to be 4/20, though 
there is no hard and fast rule that says that we must adhere to 
this last condition, which is dictated more by the numerical 
discretization of the integral forms of Maxwell’s equations 
(3), than by anything else. At this point, we note that since 
the dimensions of region-2 in the physical and virtual 
domains are different, we typically choose the radius a such 
that a << b, in order to ensure that the scattering from the 
small cylinder in the virtual domain would be small—in fact 
vanishingly small in the ideal case—in order to render it 
invisible as a—0 in the limit. Note that we must choose a 
non-uniform mesh in the virtual domain, so that we can 
simultaneously satisfy condition (i), as well as the 
constraint on the number of cells in the radial direction in the 
two domains, namely that they be equal. Though we have 
some flexibility in terms of the variation of the cell size in 
the radial direction in the virtual domain, we choose this 
variation in Ar to be smooth, and monotonically increasing in 
terms of the cell size as we go from r = a to r =c, so that the 
summation of all the Ar s equal (c —a). An example of such a 
mesh is shown in Fig. 5. 


2c’ 







Region 2: Cloak Material 


Region 1: PEC 


(a) (b) 


Fig. 5. Cloak problem with mesh schematics: (a) physical and (b) virtual 
domains. 


Having defined the meshes in the two domains, we finally 
turn to the task of determining the material parameters of the 
cloak in the physical domain. Recall that we wish to impose 
the criterion that the two sets of fields, namely (E4, H1) 
and(E-, H) in the two domains, respectively, be identical. 
Eqn. (2) tells us that the (€; u) values in the physical 
domain, must be €&/y, and uoy, where y„is the ratio of the 
areas of the n”cell in the physical and virtual domains, 
respectively. It is evident that, under these conditions, the €; 
and w;must start out at &/y,; and uo/yı, where y;the ratio of the 
dimensions of the first cell (at r = b) in the physical domain 
to that of the dimension of its counterpart (at r = a) in the 
virtual domain. Also, we recall from our previous discussion, 
that € and yw values are identical in both domains when we 
reach r = C, and that they are both just (€o, 4o), i.e., material 
parameters of free space. 


We now make two important observations. First, the cloak 
we have designed by following the procedure just described, 
which is based on the TO algorithm—though implemented 
differently than via the use of the Jacobian—varies 
inhomogeneously in the radial direction, albeit smoothly. 
Second, the ¢ and u values in the physical domain are larger, 
by a factor of //y;, than their free-space counterparts. This 
factor is responsible for the root cause of difficulty 
encountered when designing TO-based cloaks, since it calls 
for Metamaterials to fulfill the requirements on the material 
parameters if we insist that the cloak must render the target 
totally “invisible”. This is because we can fulfill that 
condition if and only if a—0, in the strictest sense, and this 
implies that y;ymust also follow suit and tend to zero as well. 


We should mention that although we chose the simple 
geometry of the cylinder to identify the fundamental 
difficulties with the implementation of the TO paradigm, the 
problem with the realization of material parameters persists 
regardless of the geometry of the target, as long as we insist 
that it becomes invisible, which in turn requires that the scale 
factor (equivalent to y) tends to 0. 


Before we discuss our strategy for overcoming this 
fundamental roadblock, we examine another fundamental 
limitation posed by the TO paradigm when we attempt to 
reduce the thickness of the cloak t (= c— b),to realistic 
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values, e.g., a small fraction of the wavelength. We note that 
the material values (£o, 4o) near the outer boundary of the 
cloak, located at r = C, do not change, regardless of whether 
r is large or small, and, for that matter neither do the 
behaviors of (¢,2) of the cloak in the neighborhood ofr = b, 
where the above parameters —œ. While we can attempt to 
partially mitigate this problem by imposing a cap on the 
values of these parameters, we cannot change the fact that 
the relative ¢, and u, must reduce from very large values at 
the surface of the cylinder to unity within a relatively short 
distance; and this, in turn, again poses realizability problems 
in practice. For this reason, practical realizability of thin 
invisibility cloaks have not met with too much success in the 
past, and it is unlikely that the thickness issue will be 
resolved anytime soon if we continue to impose the above 
invisibility criterion on the TO-based cloak designs. Fig.6 
shows the deterioration of the performance of the cloak (see 
Fig.6(c) when we design a thin three-layer cloak by using 
discretized values of the cloak parameters, which vary 
continuously from r = b to r = c). It is evident that the 
wavefront of the field becomes considerably distorted when 
we compare the performance of the thin cloak with the ideal 
one the results for which are included in Fig.6 (b) for 
comparison. 


(a) 





Fig. 6. Field distribution for A PEC cylinder in (a) an ideal TO cloak, (b) 
an ideal, thin TO cloak, (c) a 3-layered TO cloak in which the medium 
parameters at each layer correspond to those of the ideal thin cloak. 


We now proceed to outline a strategy for designing 
coatings for radar targets which circumvent the various 
problems identified above, when attempting to realize 
realistic cloak designs. 


II]. RCS REDUCTION FOR PRACTICAL TARGETS BY USING 
A MODIFICATION OF THETO CONCEPTS 


We now present a modified version of the TO approach 
that attempts to overcome the major problems encountered in 
TO-based designs. We make it clear, at the outset, our goal is 
to reduce the RCS, rather than attain invisibility. 


Our strategy is to follow a three-step process, which we 
will now outline below. Our first step is to design wideband 
layered absorbers for infinite, planar, conducting ground 
planes. We note that this is alsothe basic approach to 
designing coating for radar targets to render them stealthy, 
regardless of their shapes and that it is relatively easy to 
carry out this design by using optimization algorithms to 
determine the material properties and the layer thicknesses to 
reduce the reflection from the coated PEC plane, say, below 
-10dB level, and over a wide frequency band, say covering 
the radar frequencies from 2 to 18 GHz, as typical examples. 
(Note that this is in contrast to 2-3% bandwidth of cloaks 
designed by using the legacy TO algorithm). 


A typical cost function [40] to be minimized, which is 





suitable for optimizing the layer thicknesses while balancing 
the absorption performance and overall thickness can be 
defined as follows: 


k 
F=m*R,+(1-m)*) d: 

i=1 
(4) 
where d; represents the thickness for each layer of the 
absorber; Rz stands for the reflection level at the air-absorber 
interface and is a function of the material parameters and 
thicknesses d; of the absorbing layers (see Fig.7); and m is 
the weight of the reflection level in the optimization, 
signifying the importance of the reflection level over the 
collective thickness of the layered absorber. Note that m = / 
corresponds to the case where the total thickness is 
predetermined, and only the reflection coefficient is 
minimized. Note that the material parameters are dispersive, 
and hence the optimization to determine the thicknesses 
should be carried out over the entire desired frequency band. 


PEC 


Fig. 7. Model schematics of the optimization problem for the multilayer 
absorber. 


Besides the wide bandwidth, one of two other important 
attributes which distinguish this type of coating design from 
the traditional TO cloaks is that the proposed blankets are 
relatively insensitive to polarization, mostly due to the fact 
that we have a wide range of homogeneous, isotropic 
materials available to be selected as absorbing layers. The 
incident angle of the impinging radar signal interrogating the 
target has a much wider tolerance range than that of the 
TO-based cloak. This is in contrast to the TO-designed 
cloaks [9], which become rapidly ineffective if either the 
polarization or the incident angle deviates even slightly from 
that of the combination of the two parameters for which the 
TO-cloak was designed. 


The second attribute of the present design is perhaps even 
more crucial than the first one we mentioned above, and this 
concerns the availability of the materials needed to design 
the blanket, which, as we might recall, posed serious 
realizability problems, owing to the unrealistically large (or 
small) parameter values (see Fig.2) demanded by the TO 
algorithm. 


As is well known, the RAMs (radar absorbing materials) 
have been around for a very long time, some for many 
decades, dating back to when stealth aircrafts came into 
vogue in the sixties, and earliest theoretical and experimental 
work starting around 1930s. We realize that information on 
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some of these RAM materials is not available openly 
because of their “classified” or “secret” nature, 
understandably so because they are used in military 
applications to design stealth aircrafts and missiles. 
Nonetheless a plethora of information about similar 
absorbing materials is available in the open literature, 
including the details of their fabrication, which have been 
described in [27, 28], for instance. 


Here we will use two different types of materials namely 
CoFe Nano-Flakes (NF) and CoFe Nano-Particles (NP), 
whose frequency variations are shown in Fig.8. We point out 
that these materials can be realized with relative ease, as is 
evident from [27, 28], where the details of their fabrication 
can be found. 
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Fig. 8. Real and imaginary parts of the (a) permittivity and (b) 
permeability for the two types of absorbing materials used in the multi-layer 
absorbing blanket design. 


We note that the e and u values of these materials are both 
complex and dispersive. We hasten to point out, however, 
that unlike the MTMs, which typically have variations with 
frequency because of resonant inclusions they use to achieve 
effective and u values that are very high, very low, or even 
negative, the material parameters of the NP and NF materials 
vary relatively smoothly with frequency, and this is crucial 
for realizing wideband performance, as we will soon see. 


To illustrate the fact that we can indeed achieve wideband 
performance in terms of reflection reduction over a wide 
frequency band with relative small thicknesses of 2, 4, 6 and 
7 layer absorbers we refer to Fig.9, in which a 10dB (or 
better) reduction in the reflection coefficient is presented. 
Although not shown here, the results for the reflection 
coefficient reduction are also satisfactory when either the 
polarization, or the incident angle is varied and this is also 
true when both are changed simultaneously. 
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Fig. 9. Reflection coefficient of multilayer absorber backed by a PEC 
plate for different number of layers. 


We now move to the second step in our design procedure, 
which is to adapt the blanket designed for the infinite PEC 


plane to an arbitrarily shaped object. Initially we consider an 
object with a smooth surface whose radius of curvature is 
moderate-to-large everywhere. We will generalize the 
procedure in the third step, using the principles of the TO 
when the above assumption regarding the smoothness of the 
object is not valid, as for instance when the object has sharp 
edges or bumps, as a general target would in practice. 


When the object has a relatively smooth geometry, we 
initially wrap the multilayer absorbing blanket, which we 
have designed earlier for the planar surface around the PEC 
target whose scattering cross-section we are attempting to 
reduce, and test the effectiveness of the blanket for the new 
object. For a wide variety of targets we have examined, a 
number of which are shown in Fig.10, we have found that 
the blanket does reduce the monostatic as well as the bistatic 
radar cross-section in the “reflection” region near the surface 
of the object for different angles of incidence and 
polarizations of the incoming wave. The results for a 
two-layer absorber are shown in Figs.11 and 12 for a 
rectangular cylinder of finite length, which we have studied 
as a test case. 
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Fig. 10. Back-scattering RCS of PEC objects: (a) plate, (b) pyramid, and 
(c) cylinder covered with 2 and 7-layer absorbers. 
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Fig. 11. Phase behavior of the E-field near the rectangular PEC cylinder, 
which is wrapped around by an absorber blanket, normally incident on the 


cylinder. 
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Fig. 12. Phase behaviors of the E-field near the rectangular PEC cylinder, 
which is wrapped around by an absorber blanket, for an obliquely incident 
plane wave. 


A simple test, which is typically applied to cloak designs, 
is to examine the wavefront of the total (incident+scattered) 
field, and see how the level of distortion of the wavefront 
decreases when the scatterer is covered by the layered 
absorbing coating. We present the plots of these wavefronts 
of the total fields for normal and oblique incidence cases for 
both polarizations in Figs.11 and 12, respectively. 


We observe that the object, which is a finite cylinder of 
heightl8cm, generates distorted wavefronts owing to the 
contribution of the scattered field from the object, even when 
covered by a two-layer blanket, designed for the infinite 
planar PEC object, for the nominal frequency range of 4.6-18 
GHz, with a nominal reflection coefficient of -10dB or less. 
However, we also note from Fig.11 that the distortion in the 
phase front is relatively small once we go above the low-end 
of the design frequency range, viz., 4.6 GHz for the planar 
geometry, confirming that the planar design performs 
reasonably well even though we are dealing with a 
rectangular cylinder now. We hasten to point out that the 
results presented in Fig.11 are not for a coating which has 
been optimized for the object at hand, and we expect some 
compromise in the performance of the coating. However, we 
can improve this performance by optimizing the parameters 


of the two-layer design, specifically the relative thicknesses 
of the two-layer, even as we maintain the total thickness 
intact. We expect the changes to be relatively minor, 
however, except for the corner regions and, hence, the 
optimization process should be realistic as well as 
numerically feasible. The above remarks are also applicable 
to the oblique incidence case, for which some sample results 
are presented in Fig.12. 





Fig. 13. Target detection in conventional radar scenarios: mono-static 
and bi-static scheme. 


It is important to point out that the strategy for designing 
the absorptive coating, presented herein, is very different 
from that employed for ideal traditional TO cloak, since the 
latter is designed to render the (object+cloak) composite to 
have a zero scattering cross-section in all directions, whereas 
the blanket design introduced here seeks to reduce both 
backscattering and bi-static scattering scenarios but only in 
the reflection region, and not in the forward-scattering 
direction. We hasten to point out, however, that this type of 
performance is perfectly well suited for modern radar 
systems, symbolically depicted in Fig.13 above, where only 
the scattering in the reflection region is of concern. 


For the final step, we consider the problem of absorber 
design when a shape perturbation is introduced in an object. 
Let us say that our modified target is the same rectangular 
cylinder we just considered above, except for a bump on the 
top surface. The extra corners introduced by the perturbation, 
be they smooth or sharp, would obviously introduce 
additional distortions in the planar phase front, and 
potentially increase the scattering level. Our objective here is 
to restore the field behavior so that it is close to that of the 
original object that we had prior to the introduction of the 
perturbation. 
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(a) 


(b) 


(c) 


oi 


(1) (2) 

Fig. 14. (a) Original 2-layer absorber wrapped around a rectangular 
cylinder with shape perturbation (left); (b) 2-layer absorber with 
TO-modified material properties around regions with the shape perturbation 
(left); (c) 2-layer absorber with TO modified material properties for all 
regions (left) and 2-layer absorber wrapped around a rectangular cylinder 


(right). 


We now outline the procedure for the blanket design for 
the new object, shown on the left in Fig.14, 1.e., 14(1), which 
is in the physical domain, and is a modified version of the 
one shown in the right side of the same figure; 1.e., Fig.14(2), 
which corresponds to the virtual domain. Note that unlike the 
cylinder example we discussed earlier, the medium in the 
virtual domain, surrounding the object, is no longer 
free-space, as was the case shown in Fig.5. Note also that the 
dimensions of the objects in the two domains are 
comparable, and are totally different from the legacy 
TO-design case, in which the scale factor between the 
dimensions of the object in the physical and virtual domains 
tends to infinity to render the target invisible. 


To find the parameters of the cloak for the modified 
geometry in Fig.l4, we revisit the integral forms of 
Maxwell’s Equations, presented earlier in (2a) and (2b), to 
relate the material parameters associated with the two 
systems shown in the figure. Fig.14 (1) shows the physical 
system with locally modified medium parameters for the 
perturbed object, while Fig.14(2) depicts the virtual system 
with the original medium parameters covering the 
unperturbed object. The field distribution near the perturbed 
object would obviously be different from that of its 
unperturbed counterpart, since the perturbation introduces 
additional scattering to the incoming wave. We link the 
change in the field with the modification in the geometry, 
and then compensate it by the changing the medium 
parameters, which are illustrated by a slight change of colors 
from the ones representing the original planar design in 
Fig.14(1.a). Given the specific profile of the perturbation and 
the simulated field distributions, the only unknowns, namely 
the medium parameters, can be derived as follows: 


u B f Ecd B pf E,-dl, 
Cfo oe ee oe 
—= ff Hy: ds; — = ff Hz dS, 
$ Ez ` dh 
= > SS 
—= ff Hz aS; 
ð 
— = ff He "dS, , $ E,» ‘dl, 
. ð 
fE,- dl, -Z Jf H, + dS, 
ð 
: — G Ey dl, — 5, SJ Az: aS: n 
17 Er.) 3r, 2 ™2 
$ E, dl, —= ff Hz dS, 
(5a) 
and 
ð 
Z _GHy dl, Fae S Be dSo 
© $H'dh +2 ffe, ads * 
(5b) 


where we have used scalar quantities for simplicity in (5) as 
though the geometry under consideration is two-dimensional, 
which it is in the present example, and we would need to 
replace the field quantities with vectors and the material 
parameters with tensors for the general 3D case. 


The field behaviors for the perturbed object with locally 
modified material parameters can be seen from Figs.15 and 
16. Fig. 15 shows that the amplitude of the scattered E-field 
is reduced, and that the phase front of the total E-field is 
approximately restored as well in the reflection region.Fig.16 
compares the amplitudes of the scattered electric fields for 
five different scenarios listed in the figure. 





(a) (b) 


Fig. 15. Phase behavior of scattered E-field for (a) perturbed object 
wrapped by blanket with original medium parameter; (b) perturbed object 
wrapped by blanket with locally modified medium parameter. 
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Fig. 16. Comparison of the amplitudes of electric fields scattered by 
different objects. 


We should clarify the fact that although we are referring to 
this geometry as a slab, what we are really dealing with is a 
wide rectangular cylinder, with a small thickness. 


We note that the modified slab does introduce additional 
scattering, and that the absorber does help reduce the same. 
We also note that the modified absorber improves the 
performance over the initial one, but only slightly, which 
shows that the planar version of the cloak is not all that 
inferior to the one modified for this type of geometry. 
Additional optimization of the modified cloak is expected to 
improve the performance even further, if so desired. 





Fig. 17. Possible domain decomposition of two aircrafts: F-16 Falcon 
fighter jet (left) and Predator Drone UAV (right) for scattering-reduction 
treatment. 


For an arbitrary target, we can first decompose the 
geometry of the target in a manner illustrated in Fig.17 into 
several large “blocks” that closely resemble previously 
investigated objects, and wrap each part of the target with an 
absorber, designed by using the approach based on shape 
perturbation of a related smooth object. We then follow the 
methodology we have described above to determine the 
material parameters of the coating. It should be evident that 
this is a far more realistic approach than transforming the 
geometries of these complex objects into an infinitesimally 
small-size target, as called for by the TO algorithm for cloak 
designs, and following the TO recipe corresponding to such a 
geometry transformation, which is bound to lead to 
unrealistic and impractical material parameters and/or 
structural elements. 


Modifying the permittivity and permeability parameters of 
the lossy materials is not the only approach one could take to 
reducing the scattering from objects that have curved 
profiles, especially when laboratory resources to fine-tune 


the material parameters in accordance with (5a) and (5b) are 
not available. 


Changing the composition of the absorber can be quite 
effective for enhancing the absorbing performance of the 
coating. Fig.18 shows the RCS results of a PEC cylinder of 
radius 3cm, illuminated by a plane wave. The bi-static RCS 
evaluated at 10GHz, which corresponds to a wavelength of 
3cm. Various thickness compositions are analyzed to 
determine the one most suited for this particular cylinder, 
while the total thickness of the absorber is kept unchanged at 
3mm. 


Fig.19 shows the result for another example of a PEC 
sphere of radius 30cm illuminated by a plane wave. The 
bi-static RCS is again calculated at 10GHz. Various 
thickness compositions are analyzed to determine the one 
that provides an optimal performance. Once again, the total 
thickness of the absorber is kept unchanged at 3mm. 


We note from the figures that the optimal choice for a 
planar absorber is not necessarily the best one for a curved 
object. A moderate performance enhancement in terms of 
absorption can be achieved by simply changing the thickness 


of the absorber, as we might expect. 
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Fig. 18. Comparison of bi-static RCS for a long PEC cylinder coated 
with absorbers with various thickness compositions, for horizontally (up) 
and vertically (down) polarized illumination. 
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Fig. 19. Comparison of bi-static RCS for a PEC sphere coated with 
absorber with various thickness compositions. 


It should be pointed out that layered absorbing materials, 
which comprise of more than two types of lossy components, 
and/or more than two layers, are likely to be better suited for 
the task of coating objects with a curved profile. Fig.20 
shows the example of a five-layer absorber designed for a 
planar surface using the optimization techniques discussed 
previously. Wideband absorption performance is achieved 
from 2 to 16.9GHz with a total thickness of 4.26mm (as 
opposed to a frequency range of 4.6 to 18GHz and a total 
thickness of 3mm) [39]. The increased number of layers 
provides increased degrees of freedom, which could 
potentially achieve the required level of performance without 
having to resort to synthesis of new material. 
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Fig. 20. Calculated and measured reflection coefficient of a five-layer 
absorber. 


Optimization techniques can still be efficiently applied to 
design absorptive coatings for objects with cylindrical or 
spherical shapes, since the fields scattered from these 
canonical shapes can be analytically expressed as functions 
of the material parameters and compositions of the coatings. 
When the shapes of the objects are no longer canonical, as in 
the case in real-world scenarios, applying optimization 
technique may be a time-consuming task. This is because the 
forward problem, namely solving for the electric field 
distribution through numerical simulation—since analytical 
calculation is not an option here—needs to be carried out for 
each and every object, for which we are designing the 
absorber. It is typically computer-intensive, and costly both 


in terms of CPU time and memory. Thus the use of powerful 
optimization algorithms is highly recommended for 
designing these absorbers. 


Fig.21 show the results for a PEC ellipsoid whose major 
and minor axes are 15cm, 3cm and 3cm respectively and 
which is illuminated by a plane wave with horizontal and 
vertical polarizations. The bi-static RCS is calculated at 
10GHz, which corresponds to a wavelength of 3cm. Two 
absorber compositions are analyzed to investigate the 
performances of these two configurations. The first one of 
these uses the original absorber design for the planar 
structure, while the second corresponds to the optimal 
thickness composition we presented in Figs.18and19. The 
total thicknesses of the two configurations are kept 


unchanged and it’s 3mm in both cases. 
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Fig. 21. Comparison of bi-static RCS for a PEC ellipsoid coated with 
absorber of various thickness compositions, under horizontally (up) and 
vertically (down) polarized illumination. 


It is evident from Fig.21 that the modified composition 
exhibits a slightly enhanced performance for both 
polarizations. However, we should point out that there is no 
guarantee that an improvement can always be achieved over 
that of the original planar design, and that an optimal design 
for the absorber coating can indeed be derived, for a given 
object by simply varying the layer thicknesses using a 
trial-and-error approach. An alternative technique is to 
modify the material parameters, as illustrated above, which 
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provides a more straightforward solution, especially for 
geometries shaped like the ellipsoids. 


Reduction of RFI, another important application for the 
absorber-wrapping is to mitigate the problem of antenna 
blockage in a shared-platform environment, as shown in 
Fig.22, in which the introduction of the aggressor antenna 
can raise the far-end sidelobe levels of the parabolic dish 
significantly. To mitigate this effect, we can wrap the 
monopole by using a multilayer absorber optimized for a 
planar geometry, as we have done in the past, or we can 
optimize the layer thicknesses and material parameter for the 
circular cylinder geometry. Alternatively, we can use a 
conducting saucer-like structure wrapped around the 
monopole, as proposed in [29-31]. Fig.23 shows the model 
used to construct the alternative cloak with the monopole as 
the cloaked target. Fig.24 shows the far field patterns of a 
parabolic dish antenna—the victim antenna—when a 
monopole, which is the aggressor antenna, is placed in the 
vicinity of the dish. The patterns of the dish/monopole 
composite are also included in the figure for two treatment 
plans applied to the aggressor antenna. We note that the 
introduction of the absorber treatment reduces the sidelobe 
levels of the dish antenna, as compared to the case for the 
far-end antenna combination without the treatment. 


However, we also note that an elaborate cloak design, 
shown in Fig 23, isn’t really needed, and same level of 
performance can be achieved with a thin absorber type of 
coating. 


oN 








500mm 
770mm > E 





Fig. 22. Model schematics for the blockage problem: (left) side view and 
(right) isometric view of the victim antenna (parabolic dish) and the 
aggressor antenna (monopole) sharing the same platform. 
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Fig. 23. Model schematics for the alternative cloak. 
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Fig. 24. Radiation patterns for the antenna blockage problem: (blue) dish 
antenna, (green) dish antenna and the monopole antenna in the vicinity and 
(red) dish antenna with absorber-treated monopole antenna. 


Recent advances in material engineering and research in 
graphene-based absorbing materials also broadens the choice 
of the materials used for the absorber designs discussed in 
this paper. Graphene-based absorbing materials generate 
sufficient level of magnetic loss without the need to 
introduce magnetic elements such as Co, Fe and Ni. These 
new materials also have the advantage of being lightweight, 
and thin, that are highly desired attributes for airborne 
applications. There have been reported cases of achieving 
below -70dB reflection level shielding with 2.09mm of 
graphene nanoplatelet-synthesized material [32], as well as 
below -20dB reflection level over the entire 4.5-18GHz with 
8mm of multilayered carbon nanotube fabric [33]. 


The performance of the absorbers can also be enhanced by 
introducing FSS-type structures within the layers, which add 
virtually no extra thickness to the composite structure. The 
FSS structures can be tailored to enhance the absorption in 
target frequency bands [34] or to broaden the working 
bandwidth [33] of the initial layered absorber. 


IV. CONCLUSIONS 


In this paper we have presented a strategy for designing 
cloaks and absorbing coatings for scattering reduction which 
deviates from the classical TO algorithm. We _ have 
discussed the limitations of the Transformation Optics-based 
algorithm for designing cloaks, which typically calls for 
Metamaterials that are not available in nature and suffer from 
problems of narrow bandwidths, high losses, unacceptable 
levels of polarization and incident-angle sensitivities, etc., 
when synthesized artificially. We have provided an 
alternative approach to TO for relating the material 
parameters of the physical and virtual domains, by using the 
integral forms of Maxwell’s equations. We have used this 
new perspective to identify the root cause of the difficulty 
that arises when designing cloaks by using the TO, namely 
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inability to realize extreme values of the material parameters 
called for by the TO algorithm 


Following this we have shown how we can develop 
absorbing coatings for planar conducting structures, and then 
use it as a stepping stone for designing absorbing blankets 
for objects of arbitrary shape. 


Finally, we have detailed an algorithm for deriving the 
material parameters of absorbing blankets for arbitrarily 
shaped objects, to reduce their scattering levels in the 
reflection region, over a wide frequency band, and for wave 
incidences with arbitrary polarization and incident angles. 
An alternative approach to change the thickness composition 
of the absorbers has also been discussed. A number of 
examples have been included to illustrate the performances 
achieved by the absorbing coatings designed by using the 
developed technique. 
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